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An implementation of premixed equilibrium chemistry has been completed for the 
OVERFLOW code, a chimera capable, complex geometry flow code widely used to pre- 
dict transonic flowflelds. The implementation builds on the computational efficiency and 
geometric generality of the solver. 
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0 angle along sphere surface(degrees) 

£, dimensionless wedge axial distance (x/L) 

c isentropic sound speed 
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L wedge length 

p static pressure 

s entropy 
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2 ^ 3xj ^ 3xi ) 

T static temperature 

u, v, w cartesian velocity components 

x, y, z cartesian position coordinates 

Introduction 

Building on the successful application of the perfect 
gas version of OVERFLOW to hypersonic flowfields, 1 
the next step in implementing finite rate chemistry 
capability in OVERFLOW is the generalization of the 
gas model from single species perfect gas to allow pres- 
sure and temperature as arbitrary functions of den- 
sity and internal energy. The specific thermodynamic 
model used was that of Liu and Vinokur. 2 This model 
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provides pressure and temperature for an equilibrium 
mixture of air over a wide temperature and pressure 
range. 


Method 

The Reynolds averaged Navier-Stokes equations, 


written in conservation law form are 
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where, we assume a simple Boussinesq type stress- 
strain relationship for the turbulence, 


Ty = 2(|i.+ p. T )Sij - jpxSij (6) 


The thermodynamic relation for pressure p = p(p, e) 
is used in place of the perfect gas assumption(p = 
p(y — 1 )e) in computing the fluxes of the momentum 
and energy equations, and the thermodynamic relation 
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for temperature T = T(p, e) enters into these equations 
indirectly, in the evaluation of molecular viscosity and 
thermal conductivity. 

Implicit solution requires the calculation of Euler 
Jacobians (f^-, , and *jq-) and for diagonalized 

methods, their eigenvectors and eigenvalues. The Eu- 
ler fluxes are simply the Navier-Stokes fluxes with 
4 , 47 -, and k, set to zero. This generally involves 
an evaluation of the isentropic sound speed, c. The 
sound speed is derived from the partial derivatives of 
the pressure returned by the thermodynamics mod- 
ule, and it current implementation allows a capability 
more general than mixtures of perfect gases in chemi- 
cal equilibrium. 

This solver is actually capable of predicting flows 
of any simple compressible substance where pressure 
and temperature are functions of only density and in- 
ternal energy. The sound speed is given by the general 
relation for isentropic sound speed c 2 = |^l As we 


have a table lookup routine which returns values of p, 


T, 


<Lr 

ap 


iE 

3e 


3T 

ap 


U , we need to be able to get 


the sound speed in terms of these variables and their 
derivatives. 

From the Gibbs equation for a simple compressible 
substance, 


ds = y (de +pdv) 


= -(de-pdp/p 2 ) 


we get the relations 
0e 
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which, along with the trivial relations 
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allows us to transform the partial derivatives we have 
naturally as a result of the bicubic spline: 

into derivatives where [p, s] are the independent vari- 
ables. Thus we may calculate the sound speed from 
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This is the relation used to give the sound speed, which 


is implemented by calculating an effective y, 


y = pc 2 /p 

£ 9p + 1^P 

P 3P £ P p 


( 11 ) 


Implementation 

The base code used is OVERFLOW 3 which is a 
finite difference, chimera capable, complex geometry 
flow code widely used for perfect gas prediction. Ma- 
trix dissipation is used with central space differenc- 
ing and multigrid to provide a rapid, robust solution 
method suitable for high speed flow with strong shock 
waves. 1 

OVERFLOW was modified to replace the perfect 
gas model built into the code with an arbitrary gas 
model, where the pressure and sound speeds are func- 
tions of the density and internal energy, as given in 
Liu. 2,4 An efficient table lookup method returns the 
pressure and temperature and their derivatives corre- 
sponding to a given value of [p, e]. A bicubic spline is 
used to produce smooth interpolations of both func- 
tions and their derivatives. 

Four additional variables are added to the conserva- 
tive variables [p, pu, pv, pw, pe]. The four additional 
field variables are [p/Poo, T/Too, Pr, p/(pRooT)]. These 
variables are obtained from the conservative variables 
at the start of each iteration, before the laminar viscos- 
ity is calculated. The perfect gas assumptions are then 
replaced by the use of the calculated pressures(in the 
momentum and energy fluxes, and their derivatives), 
temperatures (in the viscosity and conductivity cal- 
culations), and sound speed (in the implicit jacobian 
calculations). 

2nd order central flux differencing combined with 
matrix dissipation 5, 6 is used in forming the right hand 
side residuals, similar to the scheme used for perfect 
gas solutions. 1 The implicit solution method used is 
the Pulliam-Chaussee diagonal algorithm, with the 
eigensystem 7 the same as used to implement the ma- 
trix dissipation. Roe averaging is used in forming the 
matrices in the dissipation scheme. 

The multigrid method 6 has been modified to en- 
hance the robustness of the method. The original 
multigrid method used a linear interpolation of the 
conservative variables to map the coarse grid solu- 
tion onto the next finer grid. This interpolation was 
replaced with a linear interpolation of the variables 
[p,u, v,w, e], which ensures a positive value of e on 
the finer grid if e is positive on the coarser grid. 


Results 

Inviscid Wedge Flow 

This flowfield is a self similar flowfield useful for 
gauging the ability to capture shocks at oblique an- 
gles to the grid, and to obtain the external inviscid 
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Fig. 1 10°M=25 Inviscid Wedge Grid, Colored by Axial Velocity 



Fig. 2 Inviscid M=25 Wedge, Profiles of Pressure, Total Enthalpy, Temperature and y 


flowfield downstream of the shock. The case consid- 
ered here is a 10° wedge, at M=25.083, with freestream 
static conditions of Poo = 20.24Pa.Too = 257° K. 

A converged solution is shown in Figure 1. The cap- 
ture of the oblique shock is clearly demonstrated, with 
uniform conditions obtained on both sides of the shock. 
The shock crosses the grid system obliquely, with the 
shock moving out one grid line for every two gridlines 
moved downstream. 

Profiles of pressure, total enthalpy, temperature, 
and pressure are shown in Figure 2. The pressure 
profile shows minimal overshoot, with good similarity 
between the various axial locations The total enthalpy 
is well conserved (for this flowfield it is constant), with 
errors, even at the shock, less than 0.5%. Total en- 
thalpy errors behind the shock are less than 0.1%. 
The post shock temperature reached in this flowfield 
is 1500°K, for which there is appreciable dissociation, 
as shown by the gamma variation from 1.4(freestream) 
to 1.3(post shock). 

Inviscid Blunt Body Flow 

This case is the inviscid flow past a blunt cylinder 
(Fig 3) over a range of freestream static pressure and 
Mach number. The pressure conditions range from 
p = lPa(h=80km) to p = 10KPa(h=16km). The 



Mach numbers range from M. «, = 8(2.4km/sec) to 
Moo = 32(9.5km/sec). The solutions were obtained 
on a 101 x 129 grid suitable for computation of a per- 
fect gas, so much of the grid lies outside the shock. 

These solutions demonstrate the ability of the ma- 
trix dissipation to capture the shock over the wide 
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a) Moo — 8 


b) Moo = 16 


c) Moo = 24 


d) Moo = 32 


Fig, 4 Inviscid Cylinder Blunt Body Residual Histories 






a) M^o — 3 


b) Moo = 16 


c) Moo = 24 


d) Moo = 32 


Fig. 5 Inviscid Cylinder Blunt Body Stagnation Line State 


range of conditions representative of earth atmospheric 
reentry. The smoothing parameters are unchanged 
over this entire range of conditions, and only the M = 
32 cases required an adjustment of the CFL(halved 
from the other conditions) to obtain convergence (Fig 

4 ). 

While all the cases show a similarity in the stag- 
nation pressure profiles, (Fig 5) the variations of the 
temperature and density are more pronounced as the 
Mach number increases. At a fixed Mach number, a 
lower ffeestream pressure implies a lower post-shock 
temperature, as the dissociation is more complete at 
lower pressure. An interesting trend is evident in the 
density profiles. As the Mach number increase from 8 


to 24, the post shock density increases, but from Mach 
24 to Mach 32, there is a decrease in the post shock 
density especially at the higher pressures. 

The preservation of total enthalpy across the shock 
is good at all Mach numbers. Along with being a 
check on the fidelity of the solution, this will be im- 
portant in obtaining good heat transfer predictions in 
viscous flowfields. The pressure profiles in all cases 
are collapsed nicely when normalized by the Newto- 
nian pressure, pu^,. 

Viscous Flat Plate 

As important as inviscid shock capturing capabili- 
ties are the abilities to predict heat transfer and skin 
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Fig. 6 Flat Plate Skin Friction Predictions 




Fig. 7 Flat Plate Heat Transfer Predictions 


friction. Solutions for a laminar flat plat at freestream 
Mach numbers of 2 and 4 were obtained on 3 grids, 
1 29 x 1 29, 65 x 65, and 33 x 33. The skin friction pre- 
dictions are all nicely in agreement with theory (Figure 
6) and the heat transfer predictions (Figure 7) are sim- 
ilarly good. These predictions were obtained with the 
same dissipation parameters used to produce the invis- 
cid wedge and cylinder solutions. Also shown are the 
skin friction predictions with the lowered dissipation 
settings, which agree quite well with the ‘standard’ 
dissipation settings. 


Viscous Sphere 

This case considers a viscous flow over a lm radius 
sphere at M,*, = 17.6, p = 57.4Pa, T*, = 360°K, as 
shown in figure 8. This case is the benchmark case 
for LAURA, (http://hefss.larc.nasa.gov). The viscous 
wall boundary condition was no slip, T w = 900°K. 
The finest grid used was 65(streamwise) by 129 (wall 
normal), with a constant wall grid spacing of 0.5p.m. 
Two coarser grids were derived by this one by choos- 
ing every other point (33 x 65), and every fourth 
point(17 x 33). 

The stagnation streamlines for this case (Figure 9) 
are given in both external flow (linear in axial distance) 
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Fig. 8 Viscous Moo = 17.6, p = 57.4Pa lm Sphere 



Fig. 9 Stagnation Profiles, 


and boundary layer (logarithmic in axial distance) co- 
ordinates. The pressure profile looks very similiar 
to the inviscid case (Figure 5), and agrees well with 
the prediction of LAURA, but here the stagnation 
enthalpy and static temperature variations are dom- 
inated by that of the boundary layer variations. The 
temperature predictions show the same boundary layer 
“lumpiness” predicted by LAURA, whereas the total 

6 



oo = 17.6, p = 57.4Pa lm Sphere 


enthalpy prediction is quite smooth over the same re- 
gion. Heat transfer predictions(Figure 10) are similiar 
to those of the LAURA code. 

High Enthalpy Nozzle Expansion 

This case highlights the ability to predict flowfields 
with extreme expansion. T10-97 is one of the First 
Europe-U.S. High Speed Flow Database test cases. 8 
This case documents to flowfield in the ONERA F4 
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Fig. 10 Wall Heating Distribution, Moo = 17.6, 
p = 57.4Pa lm Sphere 



a) Density, Overall Flowfield 



b) Reservoir Detail 

Fig. 11 T10-97 (ONERA F4) Wind Tunnel Nozzle 
Flowfield 

high enthalpy facility. The resevoir conditions are p = 
373bar, hp = 140Too, and the exhaust conditions are 
Too = 273° K, poo = 68Pa. The resevoir is charged 
with synthetic air, composed of ^Oz and ^N 2 by mass 
fractions. 

A sequence of three grids was again utilized, with 
the finest grid 513(Streamwise) x 129(WallNormal), 
and the coarser grids again created via ’bicimation’. 
The solutions for this case depend critically on the 
air thermodynanmic models, 7 8 and the current results 
(Figure 12) are consistent with the results of other 
equilibrium air codes utilized to simulate this case. 

Conclusions 

The present method shows promise for rapid predic- 
tion of complex flowfields where fully mixed chemical 
equilibrium is an appropriate approximation, and is 
another step in the overall goal of adding chemically 
reacting flow capability to OVERFLOW. 
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